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ABSTRACT 

In the previous papers in this series, we found that radiative torques can play a 
major role in the alignment of grains with the interstellar magnetic field. Since the 
radiative torques can drive the grains to suprathermal rotational speeds, in previous 
work we made the simplifying assumption that the grain principal axis of greatest 
moment of inertia is always parallel to the grain angular momentum. This enabled us 
to describe many of the features of the grain dynamics. However, this assumption fails 
when the grains enter periods of thermal rotation, which occur naturally in the radiative 
torque alignment scenario. In the present paper, we relax this assumption and explore 
the consequences for the grain dynamics. We develop a treatment to follow the grain 
dynamics including thermal fluctuations and “thermal flipping”, and show results for 
one illustrative example. By comparing with a treatment without thermal fluctuations, 
we see that inclusion of thermal fluctuations can lead to qualitative changes in the 
grain dynamics. In a future installment in this series, we will use the more complete 
dynamical treatment developed here to perform a systematic study of grain alignment 
by radiative torques. 

Subjeet headings: ISM: dust, extinction — polarization — scattering 


1. Introduction 

Polarization of starlight by the interstellar medium was discovered serendipitously (Hall 1949; 
Hall &: Mikesell 1949; Hiltner 1949a, b). It was immediately recognized that the polarization 
must be due to selective extinction by dust grains, and that this required that interstellar dust 
grains be both nonspherical and aligned. The first papers seeking to account for the observed 
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alignment appeared almost immediately (Spitzer & Schatzman 1949; Spitzer & Tukey 1949; Davis & 
Greenstein 1951; Spitzer & Tukey 1951), but a satisfactory explanation for the grain alignment has 
eluded theoretical understanding for many decades, despite numerous investigations (see Lazarian 
2002 for a recent review). 

Davis & Greenstein (1951) suggested that the alignment results from paramagnetic dissipation. 
For this mechanism to act, the grain must be paramagnetic and it must rotate through the static 
interstellar magnetic field. Paramagnetism appears likely for amorphous silicate grains, and plau¬ 
sible for hydrogenated carbonaceous grains. Indeed, given the substantial fraction of the overall 
grain mass contributed by Fe, a significant ferromagnetic or ferrimagnetic fraction would not be 
impossible. 

Davis & Greenstein assumed that the grain rotation is excited by elastic impacts with gas 
atoms. In this case, the energy in rotation about any grain axis is ~ ^fcTgas, where k is Boltzmann’s 
constant and Tgas is the gas temperature; thus, this motion is called “thermal rotation”. The grain 
rotating through a static magnetic field is analogous to a magnetic resonance experiment, in which 
a static solid is exposed to a rotating field. In both cases, the material tries to magnetize along 
the field direction, but with a lag, and energy is dissipated into heat. In the magnetic resonance 
experiment, the dissipated energy originates in the radiation field, whereas in the interstellar case, 
the grain’s rotational energy is dissipated. 

The statistical mechanics of a rigid paramagnetic grain was analyzed by Jones & Spitzer (1967) 
and Purcell & Spitzer (1971). If the gas temperature Tgas exceeds the dust temperature and 
if there were no disalignment due to random collisions with gas atoms, then the Davis-Greenstein 
mechanism would drive a grain to tend to spin about its principal axis of greatest moment of inertia 
(hereafter ai), which would tend to be aligned with the local interstellar magnetic field B. 

However, for typical conditions in the diffuse ISM, the timescale for disalignment (due to 
collisions with gas atoms) is substantially shorter than the Davis-Greenstein alignment timescale 
tdg- Jones & Spitzer (1967) showed that if grains have superparamagnetic inclusions (e.g., domains 
of pure Fe) then tdg could be reduced by several orders of magnitude. Mathis (1986) noted that, 
since larger grains would be more likely to contain one or more such inclusions, this could explain 
the observation that relatively large grains (a > O.l^m) are well aligned while smaller grains 
(a < 0.05;um) are not (see, e.g., Kim &: Martin 1995).^ 

Research on the rotational dynamics of interstellar grains discovered important effects which 
had hitherto been overlooked. Martin (1971) pointed out that a charged spinning grain has a 
magnetic moment parallel or anti-parallel to the angular velocity cj, so that the grain processes 
about B. The precession rate is fast enough that, even if the actual grain alignment mechanism did 
not involve B, the observed polarization would be either parallel or perpendicular to B. Dolginov 
& Mytrophanov (1976) showed that the Barnett effect (the tendency for a spinning body to become 


^Throughout this paper, we characterize the grain size by the radius a of a sphere of equal volume. 
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magnetized antiparallel to cj) provides a much larger magnetic moment. 

The problem of grain alignment therefore cleanly separates into two issues: (1) the alignment of 
grain angular momentum J with B, and (2) the alignment of the grain itself with its instantaneous 
angular momentum J. 

Purcell (1979) realized that the Barnett effect would lead rapidly-rotating grains to tend to spin 
about ai, its principal axis of largest moment of inertia. This is the conhguration for which a grain 
with constant angular momentum J has minimum rotational kinetic energy. If the grain starts 
in a different configuration, then generally the magnetization due to the Barnett effect changes 
direction (in a periodic manner) in grain body coordinates; the resulting paramagnetic dissipation 
(which Purcell called “Barnett dissipation” in this case) drives the grain to its state with minimum 
rotational energy.^ Purcell found that the alignment of ai with J occurs on a much shorter timescale 
than tdg or the gas-drag timescale, if the grain rotational kinetic energy kT^. The Barnett 
dissipation is necessarily accompanied by thermal fluctuations which act to disalign ai from J 
(Lazarian & Roberge 1997). Roberge & Lazarian (1999) have calculated the Davis-Greenstein 
alignment of oblate spheroidal grains with this effect included. 

Purcell (1975, 1979) also pointed out that grains are subject to systematic torques fixed in 
grain body coordinates, and that these torques can spin the grain up to suprathermal rotational 
speeds. A suprathermally-rotating grain is largely impervious to disalignment by random gas atom 
bombardment, and is thus free to undergo alignment on the Davis-Greenstein timescale. However, 
the systematic torques discussed by Purcell result from processes that occur at the grain surface 
(e.g., the formation of H 2 at special surface sites, with subsequent ejection), and thus depend 
sensitively on the details of the grain surface. Changes in the surface (e.g., due to accretion of 
atoms from the gas) can change the magnitude and direction of the systematic torque. As a result, 
the grain will sometimes be spun down, when the torque is directed opposite u, and the grain 
rotational kinetic energy can be reduced to ~ /cTgas for a period of time. Such episodes are called 
“crossovers” 

Crossovers were first studied by Spitzer & McGlynn (1979). They found that, in the absence 
of stochastic torques, the direction of J remains constant. The magnitude of J decreases and then 
increases again; during this episode the grain flips over. With stochastic torques, the scenario is 
largely the same, except that the direction of J after the crossover is not identical to its initial 
direction. Spitzer & McGlynn concluded that a grain would become completely disaligned after 
passing through a small number of crossovers. This implies that if crossovers occur frequently (on a 
timscale shorter than tdg)) then paramagnetic relaxation (without superparamagnetic inclusions) 


^Even in the absence of Barnett dissipation, viscoelastic dissipation would align ai with J; Purcell found that 
Barnett dissipation is stronger. 

®The systematic torque need not be directed exactly antiparallel to J to result in a crossover; any torque component 
perpendicular to J averages to zero during the grain rotation. 
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is ineffective at aligning grains. 

Lazarian & Draine (1997) pointed out that thermal fluctuations in a grain can play an im¬ 
portant role in crossovers. In the process of Barnett dissipation, rotational energy is transferred 
to the lattice vibrational modes, which act as a thermal reservoir. Energy can also move in the 
opposite direction; a Barnett fluctuation is the spontaneous transfer of some amount of energy 
from the vibrational modes to the grain rotation (at constant J). Since the rotational energy is a 
function of the angle 7 between ai and J,^ thermal fluctuations ensure that ai is never perfectly 
aligned with J, even when the grain is rotating suprathermally. Lazarian & Draine found that 
the resulting small amount of disalignment during periods of suprathermal rotation decreases the 
degree of disalignment that occurs during a crossover. Essentially, the initial disalignment limits 
the minimum value that J assumes during the crossover, and thus limits the effect that random 
gas atom impacts can have on a grain. Lazarian & Draine concluded that ordinary paramagnetic 
dissipation can indeed align interstellar grains. 

Lazarian & Draine (1999a) identified another potentially important consequence of Barnett 
fluctuations, which they called “thermal flipping”. Barnett fluctuations cause the grain to undergo 
a random walk in the angle 7 , and if the fluctuations are strong enough then the grain can flip.® 
Once a grain that is entering a crossover has flipped, the component of the systematic torque along 
J is parallel to J, rather than anti-parallel to J, and the grain is now spun up rather than down. If 
a thermal flip occurs before J becomes very small, then disalignment can largely be avoided. In the 
Spitzer & McGlynn (1979) crossover scenario, J goes to zero during the grain’s flip; in the Lazarian 
& Draine (1999a) scenario, the grain can flip at an earlier stage and bypass the disaligning, low-J 
stage.® 

Lazarian & Draine (1999a) pointed out that if a grain can thermally flip once, then perhaps 
it can do so twice. Indeed, if the Barnett fluctuations are strong enough, then the grain could 
rapidly flip back and forth. If the systematic torque (e.g., due to H 2 formation) is fixed in grain 
body coordinates, then the torque reverses in inertial coordinates each time the grain flips. With 
rapid flipping, the systematic torque will time-average to zero, and the grain is prevented from 
being spun up to suprathermal rotation! Lazarian & Draine refer to this condition as “thermal 
trapping”. They found that grains smaller than a critical size Uc ~ O.OI-O.I/xm become thermally 
trapped, which could explain the observation that relatively small grains are not aligned. 

Thus far, all paramagnetic effects had been assumed to be due to electrons. Lazarian & Draine 


^For an axisymmetric grain with specified angular momentum and rotational energy, the angle 7 is constant. The 
torque-free motion of a grain with arbitrary shape is more complicated; see §2.5. 

®For an axisymmetric grain, a flip occurs when ( 7 r /2 — 7 ) changes sign. See §2.5.2 for a description of flips in the 
general case of arbitary grain shape. 

®Of course, a flip can also occur as a result of the random torques due to gas atom impacts. However, the timescale 
for such flips to occur is much longer than the timescale for flips due to Barnett fluctuations. Furthermore, such flips 
only occur during a disaligning, low-J stage. 
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(1999b) found that for relatively low rotational speeds (as encountered during crossovers), nuclear 
paramagnetism can be much more important in Barnett dissipation and fluctuations. Taking this 
into account, Lazarian & Draine (1999b) revised their estimate of the critical trapping grain size 
to ttc > 1/um. Grains smaller than ~ l|Um are then expected to not be aligned, in conflict with the 
observation that grains with a > 0.1//m are well aligned. 

A resolution of this problem can be sought in another line of development in grain alignment 
theory. Harwit (1970a,b) pointed out that absorption and emission of photons could lead to a 
random walk in the grain angular momentum, with a tendency for the grain angular momentum 
to be parallel to the Galactic plane. Dolginov (1972) observed that individual grains might have 
chirality, with different absorption and scattering cross sections for left- and right-handed circularly 
polarized light, so that anisotropic starlight might exert a systematic torque on a grain. 

In the first paper of this series (Draine & Weingartner 1996, hereafter Paper I), we evaluated 
the radiative torque (due to the absorption and scattering of starlight) exerted on an irregularly 
shaped grain with the optical properties of “astronomical silicate” (Draine & Lee 1984). We found 
that if the interstellar radiation field is anisotropic at a level of ~ 10%,'^ then radiative torques can 
dominate H 2 formation torques for grains with a > 0.1/um. Radiative torques are not effective at 
spinning up a < 0.05/im grains in regions where H is neutral; the strong opacity of the gas beyond 
the Lyman limit denies the small grains short-wavelength radiation with which they can strongly 
couple. 

Radiative torques have clear advantages over H 2 formation torques as an agent for driving 
long-lived suprathermal spin-up in the Purcell alignment scenario. First, radiative torques depend 
on the global grain geometry and the starlight anisotropy direction, which are expected to be 
stable for > 10^ yr, longer than the timescale for grain resurfacing. Thus, spindowns and crossovers 
might occur less frequently if radiative torques are responsible for suprathermal rotation. Still, 
even one crossover can be deadly if thermal trapping occurs. Since radiative torques depend on the 
starlight anisotropy direction, which is fixed in space rather than in grain body coordinates, the 
radiative torque vector is not fixed in body coordinates, and therefore does not exactly reverse in 
inertial coordinates each time the grain flips. Thus, grains spun up by radiative torques might be 
impervious to thermal trapping. 

Radiative torques have even greater potential; in the second paper of this series (Draine & 
Weingartner 1997, hereafter Paper II), we showed that they can directly align grains. In fact, 
we found that direct alignment by radiative torques likely dominates paramagnetic dissipation for 
grains with a > O.l^m. In order to simplify the equations describing the grain dynamics, we 
assumed in Paper II that (1) ai is always parallel to J. We also assumed that (2) the timescale 
for J to process about B is much shorter than the other timescales involved in the dynamics, so 


^Weingartner & Draine (2001) have estimated that the visible/UV radiation in the solar neighborhood is ~ 10% 
anisotropic. 
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that we could immediately average over the precession in the dynamical equations. Given these 
assumptions, the dynamical state of the grain can be completely specified by the angular speed 
LV and the angle between ai and B. We constructed “trajectory maps” for three irregularly 
shaped grains (with a = 0.2/im and astronomical silicate composition) showing the evolution of 
grains initially characterized by arbitrary combinations of The trajectories often terminate 

on stable stationary points characterized by a particular i.e., the radiative torques often 

result in grain alignment. The trajectories also often pass through cj = 0, suggesting that crossovers 
(of a sort) occur naturally in the radiative torque alignment scenario. 

However, because of the simplifying assumptions adopted in Paper II, we were unable to follow 
the dynamics through these low-J stages. First, we ignored the stochastic nature of the torques 
due to H 2 formation and collisions with gas atoms. This stochasticity can be safely ignored during 
periods of suprathermal rotation but should be included when a; —> 0. Second, the assumption 
that ai II J breaks down when J ^ 0. We would expect this simplification to be more severe, since 
thermal flipping might prevent uj from becoming so small that the first simplification fails badly. 

Here we generalize the treatment of Paper H by relaxing the assumption that ai || J. In future 
papers, we will add a treatment of the stochastic nature of gas atom impacts and H 2 formation and 
we will use the resulting formalism to conduct a systematic study of grain alignment by radiative 
torques. The outline of the present paper is as follows; 

In §2, we describe the general dynamics of an asymmetric grain and how it responds to external 
torques. Most studies of grain alignment adopt axisymmetric grain shapes for simplicity, but in 
this case the high degree of symmetry suppresses the radiative torques. Thus, we treat the most 
general grain shape, with no degeneracy in the eigenvalues of the inertia tensor. We explain the 
concept of “flip state” for the grain, and show that the grain can change flip states only when 
the rotational kinetic energy passes through a critical value. Because the tumbling motion of the 
grain is very rapid, it is necessary and appropriate to average over the torque-free motion before 
considering the effects of any external torques. The torque-free motion includes exchange of energy 
between rotational and vibrational modes through the phenomenon of “Barnett dissipation”. 

In §3, we discuss the various external torques that act on interstellar grains, due to the Barnett 
magnetic moment, gas and IR emission drag, H 2 formation, paramagnetic dissipation, and starlight. 

The rapid exchange of energy between rotational and vibrational modes implies that even when 
the grain has a fixed angular momentum J, it is important to average over the different possible 
values of the grain kinetic energy; this thermal averaging is discussed in §4. In §5, we obtain the 
equations for the time evolution of the angular momentum J and the angle ^ between J and B. 

As already noted, the grain dynamics involves the possibility that thermal fluctuations will 
cause the grain to change from one “flip state” to the other. In §6, we develop an algorithm for 
evolving the grain dynamics, including stochastic flipping. We present an algorithm that can be ap¬ 
plied both in the limit where flipping occurs many times even for the shortest feasible computational 
time step, as well as in the limit where flipping occurs very rarely. 
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In §7, we present results of the dynamical evolution for one particular case (i.e., for a given 
grain shape, composition, and size; angle between the magnetic field and the radiation anisotropy 
direction; and radiation field spectrum). For this case, the analysis of Paper II yields two stable 
stationary points and two crossover points. Applying the more complete analysis developed here, 
we find a new stable stationary point with a thermal rotational speed and no crossovers. All of 
the trajectories terminate on either the new stationary point or on one of the two stable stationary 
points with suprathermal rotation. 

In § 8 , we summarize the results of this paper and discuss the remaining issues that must be 
addressed to complete our understanding of grain alignment by radiative torques. 

Since we employ a large number of physical quantities, we provide a glossary of notation 
(Appendix E) for easy reference. 


2. Grain Dynamics 
2.1. Grain Geometry 

We consider the irregular grain shape from Paper I (“shape 1” in Paper II). The moment 
of inertia tensor has eigenvalues Ii > I 2 ^ Is, with principal axes ai, a 2 , and £ 13 . We define 
dimensionless parameters aj by 

2 n 

Ij = aj-pVa^ , ( 1 ) 

where p is the density and V = 47 ra ^/3 is the grain volume. 


2.2. Coordinate Systems 

Since we are interested in grain alignment with the magnetic field B, we adopt a coordinate 
system x-q, ysi z-q, with B || zb. This system, which we call “alignment” coordinates, defines an 
inertial reference frame. It will also be convenient to define “angular momentum” coordinates xj, 
yj, zj, with zj II J. We adopt the following transformation between these coordinate systems: 

xb = cos cos (/> xj — sin (() yj + sin ^ cos 0 z j , ( 2 ) 

yB = cos ^ sin (() xj + cos (() yj + sin ^ sin (/) zj , (3) 

Zb =-sin^xj Tcos^zj , (4) 

where ^ and (/> are, respectively, the polar and azimuthal angles of J in alignment coordinates. 
Thus, we have the correspondence xj = 4, yj = </>, zj = J. See Figure 1. 

For computing radiative torques, we define “scattering” coordinates ei, 62 , 63 , with ei parallel 
to the radiation propagation direction. Three angles are required to specify the grain orientation 



in this frame. Again, see Figure 1. The orientation of ai is described by the two angles 0 G [0, vr] 
and € [0,27r], where 


ai = cos 0 ei + sin 0 cos e 2 + sin 0 sin 03 . (5) 

A third angle (3 G [0,27r] describes rotation of a 2 about ai: 
a 2 = — sin 0 cos /? ei + (cos 0 cos cos (3 — sin d> sin (3) 62 + (cos 0 sin cos f3 + cos sin (3) 03 . ( 6 ) 
Given a* and e*, the angles 0, <1>, and (3 can be found as follows: 

0 = cos“^(ai • ei) , (7) 

_i f sin 0 — ai • 02 


= 2 tan 


and 


(3 = 2 tan 


-1 


ai • 03 

sin 0 + a 2 • Oi 


( 8 ) 


(9) 


sin 0 (a 2 • 03 cos < 1 > — a 2 • 02 sin <h) 

In Appendix A, we give expressions for and (3 when the denominators in equations ( 8 ) and (9) 


are zero. 


2.3. Thermal Rotation 


In the absence of (a) internal dissipation and (b) any external torques besides those due to 
elastic collisions with gas atoms, the rms rotation rate for a sphere of radius a is given by 


f ISfcTg;; 

WT = I - 


Sirpa^ 


1/2 


-3\ 1/2 


= 1,66 X 10= ) 


100 K 




( 10 ) 


2.4. Dynamical Equations 

The equations describing the evolution of J and the rotational kinetic energy E with time are 


and 



(m 

dt 



du> ^ 


( 11 ) 


( 12 ) 


where T is the net torque on the grain. Barnett dissipation and fluctuations are manifested in the 
J • du/dt term. It is convenient to introduce the following dimensionless quantity; 


q = 2hEl.f 


( 13 ) 
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Since the rotational energy lies between J^/ 2 /i and J^/ 2 / 3 , 1 < 9 < / 1 // 3 . The dynamical equation 
for q is 

dq 1 


/i J 


du 

dt 


+ r • a; — 2qJ ■ F 


(14) 


Generally, u) is not parallel to J, in which case u executes some sort of periodic motion in 
grain body coordinates and the grain executes a more complicated quasi-periodic motion in space 
(see §§2.5.1 and 2.5.3). The timescale for this motion is usually much shorter than the timescale 
on which J changes due to the external torques (also is usually much shorter than this latter 
timescale). Thus, we will average equations (11) and (14) over one “cycle” of torque-free motion. 

We expect the terms F • u and —2q3 • F//i in equation (14) to be of the same order of 
magnitude, and to be small compared with the rate at which the rotational energy varies due to 
internal processes, so that only internal relaxation need be considered in determining the evolution 
of q. (We use the term “relaxation” to refer to both dissipation and fluctuations.) We will examine 
the validity of this approximation on a torque-by-torque basis in §3. 


2.5. Torque-free Motion 

2.5.1. Euler Equations 


The motion of an asymmetric rigid body, when the net external torque is zero, is nicely de¬ 
scribed in §37 of Landau & Lifshitz (1976). First, the Euler equations are solved for the components 
of the angular velocity u in grain body coordinates. The solution depends on the magnitude of 
the angular momentum J and the rotational kinetic energy E, and involves the Jacobi elliptic 
functions. There are several different notational conventions for these functions in the literature; 
thus, we begin by defining ours. The elliptic integral of the first kind is given by 


F(e, m) 



— m sm d) 


(15) 


The Jacobi amplitude am(rc, m) is defined as the inverse of the above integral; if rc = F{e, m), then 
e = am{w,m). The Jacobi elliptic functions are given by sn{w,m) = sine, cn{w,m) = cose, and 
dn(rt;, m) = (1 — ms\r? e)^/^. 


Next, we introduce the dimensionless quantity 

.2 ^ {h - h){(l - 1 ) 

“ (/l-/2)(l-/3g//l) 


(16) 


When q < Ii/ 12 , < 1 and the solution of the Euler equations is as follows; 


.J 

U)l = ±-— 

li 


h - hq 


-|l /2 


h-h 


dn(r, k^) 


(17) 
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J 

UJ2 — — — 

, , J 

W3 = (±)i — 

^3 




h-h _ 
h{q-l)^^/^ 


cn(r,A:^) , 


h-h 

where the (±)i in equation (19) stands for the same sign chosen in equation (17) and 

'{h-h){l-hq/hh^/^ 


The period in time t is 


T = tJ 


Pt = 4F{TT/2,k^) 


hhh 

hhh 


{h-h){l-hq/h) 


1/2 


J 


-1 


and the period in the variable r is 


= 4F(7r/2,A:2) . 


(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


When q > h/h, 


with 


and 


CUi = i 


J_ 

h 


h - hq 
h-h 


1/2 


cn(r, k ) , 


102 — — 


J_ 

T 2 


l2{l-hqlh)V'^ 


^3 = 


T = tJ 


h-h 


sn(r,/c ^) 


dn(r, k 




h-h 

(/2-/3)(g-l)^'/' 


hhh 


Pt = AF{7rl2,k-^) 


hhh 


1/2 


LU2-/3)(g-l)J 
Pr = 4F{TTl2,k-^) , 


J 


-1 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


where (±)i in equation (25) stands for the sign choice in equation (23). Of course, when q = h/h, 
wi = W 3 = 0 and UJ 2 = PJ/h- 


In Figure 2, we plot LOi (with plus signs in eqs. 17, 19, 23, and 25) over two periods for the 
case with Ji : /2 ; 13 = 3 : 2 : 1 and several values of q. 


Note that when q h/h from below [above], equations (17) through (19) [(23) through (25)] 
do not tend towards steady rotation about a 2 . Rather, they tend towards a state that oscillates 
between steady rotation with L 02 = +J/h and 102 = —J/h, with abrupt transitions between these 
two states of rotation. This is illustrated in Figure 3, where we plot uji for q slightly less than h/h- 
In Figure 3, “+” denotes the case with plus signs in equations (17) and (19) (the positive flip state 
with respect to ai; see §2.5.2) and ” denotes the case with minus signs in equations (17) and 
(19) (the negative flip state). Note also that the period Pr ^ 00 as q ^ h/h- 
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2.5.2. Grain Flips 

From equations (17) and (23), we see that the grain state is not completely specified by (J,(?). 
Given J, there are two possible configurations for which a grain can spin about one of its principal 
axes aj—either || J and > 0 or a* || — J and u)i < 0. We will call the former configuration the 
“positive flip state” with respect to axis a* and the latter configuration the “negative flip state” 
with respect to a,. The flip states can be generalized in a natural way to cases for which J does 
not lie along a principal axis. When q < I1/I2, solutions of the Euler equations with a plus (minus) 
sign in equation (17) have wi > 0 {oji < 0), corresponding to a positive (negative) flip state with 
respect to axis ai. (Also, 102 and cus average to zero.) When q > h/h^ the plus (minus) sign in 
equation (23) corresponds to a positive (negative) flip state with respect to axis 

From equations (17), (19), (23), and (25), we see that a grain can evolve smoothly from one 
flip state to another only when wi = 0 and ws = 0. This only occurs when q passes through 
Ii/l2- Figure 3 shows w* for the two flip states about ai when q is slightly less than / 1 // 2 . When 
t/Pt = 0.25, 0.75, 1.25, and 1.75, the coi are nearly identical for the two flip states and the motion is 
very similar to instantaneous rotation about a 2 . Since u)i does not quite reach zero when q < I1/I2, 
these similarities are not quite exact. 

In Figure 4 we plot liiOi/J, both for q slightly less than I 1 /I 2 (solid curve) and for q slightly 
greater than I 1 /I 2 (dashed curve).® For these curves we have adopted Ii ; /2 : ^3 = 3 : 2 : 1 and 
positive flip states. Suppose the grain has g < / 1 /I 2 at t = 0 and evolves along the solid curve 
to q = I 1 /I 2 when t = Pt/^. At this point, several different outcomes are possible. If q remains 
exactly equal to h/h, then ui and ^3 remain zero, and 102 = — J/ 12 - If <? continues to increase, so 
that q > I 1 /I 2 , then there are two possibilities; 1 . The oji emerge along the dashed curves rather 
than along the solid curves, so that the grain is now in the positive flip state with respect to a 3 . 2 . 
uji and W 3 emerge along the reflection of the dashed curves through the lines uji = 0 , so that the 
grain is now in the negative flip state with respect to ^ 3 . 

Finally, suppose q reaches a maximum value of I1/I2 at t = Pt /4 and then decreases. Again 
there are two possibilities: 1. The oJi emerge along the solid curves, so that the grain remains in 
the positive flip state with respect to ai. 2 . oji and emerge along the reflection of the solid 
curves through the lines LVi = 0. This is a grain flip; the grain is now in the negative flip state with 
respect to ai (and the angle 7 between J and ai now exceeds 7 r/ 2 ). However, such a flip will occur 
with infinitesimal probability, since there is infinitesimal probability that the maximum value of q 
attained during a fluctuation will be exactly Ii/l2- Real flips will occur in a two-step process, in 
which the grain configuration evolves across q = I1/I2 and then back again. 

Recall that as q h/h, Pt 00 . Thus, as q gets very close to h/h, the evolution in q 
occurs essentially instantaneously compared with the timescale on which the w* change for a given 


®Note that the curve for 3 is the same for both cases. 
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q. The evolution to g = I 1 /I 2 can occur at any time. As q changes, the curves for the uji change, 
with wi = 0 and cus = 0 when q = Ii/I‘ 2 - In the examples in this section, we have chosen to have 
the evolution to g = 0 occur when cus = 0 and iOi is nearly zero for the initial q (when t = Ti/4). 
We made this choice so that only two curves would be needed in Figure 4. 

In all cases, when q crosses q = I 1 /I 2 , the grain is equally likely to emerge in either flip state. 


2.5.3. Motion of Grain Axes 

Next, we need to hnd the motion of the grain axes in space. The orientation of the grain axes 
in angular momentum coordinates can be expressed using Eulerian angles (a, 7 , C)- Suppose the 
grain axes a 2 , as, ai are originally lined up with xj, yj, zj. Any grain orientation can be obtained 
by performing the following operations on the grain coordinates: 1. Rotate through angle C, about 
ai = zj. 2. Rotate through angle 7 about a 2 . 3. Rotate through angle a about ai. (See Fig. 1.) 

By resolving J onto grain axes, we find the following relations for the Eulerian angles a and 7 ; 


sinasiny =/ 2 ‘^ 2 /<^ ; 

(29) 

cos a sin 7 = I 3 UJ 3 / J , 

(30) 

cos 7 = J . 

(31) 


Angles a and 7 are periodic functions of time, with period Pt (eqs. 21 and 27). The remaining 
Eulerian angle, can be expressed as a sum of two periodic terms, one with period Pt, and the 
other with an incommensurate period. Consequently, the overall torque-free motion is not periodic. 
However, given a function A{a, (jj), we can still find its average value A(q, ±) over the torque-free 
motion, as follows: 

= dCA{a,C,3) ■ (32) 

Pt Jo 27r Jq 

In equation (32), we have explicitly indicated the dependence of A{q, ±) on q and on the flip state; 

(“—”) indicates the positive (negative) flip state with respect to ai or ^ 3 , depending on q. Of 
course, A may also depend on other variables (e.g., J), depending on which quantity A is being 
averaged. 


2 . 5 . 4 . Barnett Dissipation 

A paramagnetic grain in steady rotation (constant u) acquires a magnetization xo^/lz 
Barnett effect), where xo is the static magnetic susceptibility and 7 g is the gyromagnetic ratio of the 
microscopic magnetic dipoles that are responsible for the grain’s paramagnetism (jg = —getiB/ft ~ 
— 1.76 X I0^s“^G“^ for electrons and 7 g = gn^N/ft ~ 1-3 x I0‘^s“^G“^ for nuclei). Thus, in 
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describing the paramagnetic response of a grain, the vector Bbe ^ ^hg may be regarded as an 
effective magnetic field, termed the “Barnett equivalent field” by Purcell (1979). 

Purcell noted that when a paramagnetic grain undergoes non-steady rotation, Bbe is not static 
in grain-body coordinates; consequently, rotational kinetic energy is dissipated. Purcell assumed 
that this dissipation rate is the same as the paramagnetic dissipation rate for a static grain exposed 
to true magnetic fields given by cj/yg, where <jJ is evaluated in the (non-inertial!) reference frame in 
which the grain is stationary. Purcell considered an axisymmetric grain and ignored the constant 
component of <^/^g lying along the symmetry axis. Lazarian & Draine (2000) showed that this 
component cannot be neglected without introducing serious error in the dissipation rate when the 
grain rotation is very fast. The Purcell prescription (sometimes modified as per Lazarian & Draine) 
for obtaining the rate at which rotational energy is dissipated is widely used in the grain alignment 
literature, but it has never actually been proved to be correct. Although errors are possible in some 
cases, we expect the prescription to be accurate when applied to the grains that polarize starlight. 

Since the angular momentum in the microscopic dipoles is generally a tiny fraction of the total 
angular momentum of the grain, we make the idealization that Barnett dissipation does not affect 
the angular momentum associated with the grain rotation; i.e., it affects the value of q but not J. 
The role of the Barnett effect is to provide a coupling allowing energy exchange (at fixed angular 
momentum) between the rotational and vibrational degrees of freedom of the lattice. 

The Barnett equivalent field for a non-axisymmetric grain is quite complicated. For example, 
suppose q < Ii/l 2 - Equations (17) through (19) reveal that, in this case; (a) The biasing field along 
ai is not constant; it has a non-sinusoidally oscillating component, (b) The rotating component 
does not consist of a single field with constant magnitude executing sinusoidal rotation. A correct 
description of the paramagnetic dissipation rate for this case would be prohibitively complicated. 

However, a detailed analysis of the Barnett dissipation rate would not be warranted, because; 
(a) The paramagnetic properties of interstellar grains are poorly known, so we cannot estimate the 
dissipation rate with high accuracy, (b) As we will see, the Barnett dissipation timescale is many 
orders of magnitude shorter than the timescales for some of the other processes that play important 
roles in the grain dynamics. Thus, we will not strictly follow the evolution of q due to Barnett 
dissipation and fluctuations, but will instead adopt a simple approximate approach, as described 
below (§6). In this approach, we will simply require an estimate of the Barnett dissipation timescale 

'T3ar- 

In estimating TBar; we will consider the case of an axisymmetric (specifically, oblate; i.e., 
Ii > I2 = I3) grain. We will derive an order of magnitude estimate of TBar and assume that 
this estimate also holds for non-axisymmetric grains. For oblate grains, Bbe consists of a static 
component plus a component oJ±l'^g that rotates in the ai — sl 2 plane, with angular 

velocity curotai- Setting I2 = I3 in equations (17) through (21), we find 

J 

uj\\ = — cos 7 , 

h 


(33) 
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and 


J 

w_L = —sm 7 , 

h 


^rot 


J{h - h) 
hh 


COS 7 


where 7 is the (constant) angle between J and ai. The energy dissipation rate is given by 


(34) 

(35) 


— “^X^^-®BE,rot^rot , (36) 

where V is the grain volume, x" is the imaginary component of the magnetic susceptibility, and 
-BBE.rot = is the rotating component of the Barnett equivalent field. 

The susceptibility can be roughly approximated using the solution of the modified Bloch equa¬ 
tions derived by Wangsness (1956, his eq. A13); 


_ X 0 '^rot ^2 _ 

1 + (wil -|- Ct;rot)^2^2 A ^‘\TiT2 


(37) 


where Ti is the spin-lattice relaxation time and T 2 is the spin-spin relaxation time. For electron 
paramagnetism, we take 70^2 ~ 10“^^(15 K/T^) s, T 2 ~ 3 x 10“^^-3 x 10“®s, and Ti ~ 10“® s 
(Draine 1996). For nuclear paramagnetism, we take xo ~ 4 X 10 ^^(15K/r), T 2 ~ 3 X 10 
3 X 10“^ s, and Ti ~ T 2 (Lazarian & Draine 1999b). 


For an oblate grain, 


1 ^ A -/3 . 2 
g = 1 H---sm 7 

h 


(38) 


From equations (33) through (38), we find that (for J held constant, as we assume to be the case 
for Barnett dissipation) 


dt ’ {l-h/h) ’ 

_-i _ 21^70^2 J‘^{h - h) _ 2’F7oT2 (Ii - h) 2 

^Bar — ^ ^ 


llD hli 


llhD h 


where uP' = J‘^ j (/ 1 / 3 ) and 


(39) 

(40) 




D = 1 -I- ^(cos^ 1 T 1 T 2 ) ~ 1 -|- ^u?TiT2 

h ^^3 


where we have taken sin^ 7 ~ cos^ 7 ~ 1/2 . Thus, 

_i 2 D 70 T 2 (Ii - h) 


to 


72/3 /3 1 + {h/ 2 h)oj^TiT 2 


(41) 


(42) 


The Barnett energy dissipation rate {dE/dt)-Ba,r r\j -{jy2h)r^^^{q - 1)(1 - qh/h)/{l - h/h). 
Note that rear is the exponential decay time for (7 — 1 when (7 — 1 <C 1. 
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Adding the dissipation rates resulting from both electron and nuclear paramagnetism, we 
estimate 

2 / 1 rr \ /o -T— 3 


^Ba’r ~ l.Tyr 


.l/im\' 


a J \ Td 


15K\ /3gem 


4.7 yr 


-1 ( UJ 


\wt 

1 


0 . 


l/rm\ ^ /3gcm 


+ 


a J 
7 X 10^ 1 


P 

-3\2 


1 


+ 


7 X 10^ 


h-h 

Is ) [l + 3xl0-6a;| ' 1 + 100w| 
15K\ / Tgas \ (h-h 

Td J VlOOKy* V ^3 


1 + 8 X 10 1 + 280p 

where ws = w/lO^ s“^ and 


P = 


u;\^/3gcm / 


UJj' 


^gas 


p ) VlOOK 


0.1pm 


(43) 


(44) 


Note the very steep dependence of the dissipation rate on a when the grain rotates thermally (recall 
eq. 10 for the thermal rotation speed wt)- For large a (for which p is small), we expect dissipation 
due to nuclear paramagnetism to dominate that due to electron paramagnetism, as pointed out by 
Lazarian & Draine (1999b). In Figure 5, we plot rear versus a for a; = cut and assuming h ~ 2 / 3 . 


3. Torques 

3.1. Coordinate Transformation 

In order to find values of r(g, ±) and F • u{q, ±) for the various torques acting on a grain, we 
need to transform vectors from grain body coordinates to angular momentum coordinates. The 
transformation is a sequence of three rotations, and is given by the following equations: 

xj = (cos a cos C — sin a sin cos 7) a2 — (sin a cos C + cos a sin cos 7) % + sin C, sin 7 ai , ( 45 ) 

yj = (cos a sin ^ + sin a cos C cos 7) a2 + (cos a cos C cos 7 — sin a sin C) % — cos C sin 7 ai , ( 46 ) 
and 

zj = sinQ;sin7a2 + cosasin7a3 + cos7ai . ( 47 ) 


3.2. Barnett Torque 

The Barnett torque Tb = MBar ^ For a grain in steady rotation, the Barnett magnetic 
moment is given by 

XoV 

hBav = -^ ; (48) 

7g 

for non-steady rotation, the magnetization is more complicated, as discussed in §2.5.4. If we assume 
that equation (48) holds even for the case of non-steady rotation, then f b can be found very easily. 
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Since u) x B = u X B, Tb = {xoV/'jgju X B. Transforming u from grain body coordinates to 
angular momentum coordinates and averaging over we find that 

u{q) = q3/Ii ; (49) 


note that cj does not depend on the flip state. Thus, 

fB{q) = ^qJ^Bsm^ , (50) 


where 


= — 


XoVB 

7g^i 


25 yr 


-1 


0 .1^m\ ^ 
a j 


3 gem 

OilP 


-3 


Xo 


3.3 X 10 


-4 


—1 

hpGJ 


(51) 


Since Hb oc xq; electron paramagnetism dominates nuclear paramagnetism here. A correct treat¬ 
ment of the magnetization for the case of non-steady grain rotation does not significantly modify the 
above result. For example, for an oblate axisymmetric grain, equation (50) need only be modified 
by the substitution® 


q^ q + {q-l)x 


l + x(2| + |-l 


1 + 



hjq-l) Ti 
h—hq T2 


(52) 


where 


X = 


( h-hq \ 

\h-h) 


1/2 



(53) 


The correction term can be safely ignored since x only becomes significant when the grain rotation 
is highly suprathermal, in which case g — 1 is small. Note that the timescale for precession about B 
is rapid compared with the timescales for all other relevant processes, except for the grain rotation 
itself and possibly the Barnett dissipation timescale. 


Since Fb-u; = (/XBar^B)-a; = (u;x/.tBar)'B, we need to consider the component of /^Bar that lags 
a; in order to find Fb • u;. The relevant component is /XBar ~ hg\ thus, Fb • u; ~ jxg- 

This is identical to {dE/dt)Bar except that the interstellar field B is substituted for the Barnett 
equivalent field to/xg- Thus, 


Fb • uj ~ 


7g,e-B 

UJ 



Xg, micB 
LO 



Bar, nuc 


(54) 


where the subscripts “e“” and “nuc” denote electron and nuclear paramagnetism, respectively. At 
the low rotational speeds that may be encountered for large grains during crossovers, the electron 
term in equation (54) may be comparable to the electron Barnett dissipation rate, but in this case 
nuclear Barnett dissipation dominates Fb • 


®We used the susceptibilities given by Wangsness (1956, his eqs. All through A13) in deriving the substitution 
(52). 
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3.3. Drag 


Collisions with gas atoms will exert a drag torque, which can be written 


Tdragjgas ~ J/7drag, 


gas 


(55) 


where, for a sphere, 


^drag,gas 


8.72 X 10^ yr 


3 gem 


-3 


— ) 

0.1//m J 


f Tgas \ / 3000cm-3K \ 

\100Ky \ nnTgas / 


(56) 


For nonspherical shapes, the drag time will tend to be somewhat larger than given by eq. (56), but 
only by factors of order ~ 1.5. The drag torque will also depend on the orientation of the grain 
relative to J, but this dependence is quite weak and will be neglected. 


There is also a drag torque due to the emission of infrared photons by the grain. This torque 
has the same form as the gas drag torque but is characterized by a drag timescale 


Tirag,em — 1-1 ^ 10 yr (Qabs) 


-1 


P 


3 gem 


-3 


O.l/xm 


'Td \ / ^ISRF \ 

15 K y y ^^rad / 


(57) 


(Paper I). In equation (57), Urad is the energy density of the incident radiation field (that heats the 
grain) and uisrf = 8.64 x 10“^^ergcm“^ is the energy density in the average interstellar radiation 
field (ISRF) in the solar neighborhood, as determined by Mezger, Mathis, & Panagia (1982) and 
Mathis, Mezger, & Panagia (1983). The grain absorption cross section is 7ra^(3abs(«) ^) and the 
angle brackets denote averaging over the incident spectrum. 

The total drag timescale is given by 


hirag 




.-1 


— 7“ -I- 7“ 

drag, gas ‘ drag, em 


(58) 


Also, 

Pdrag ■ ^ — 7i(U^/Tdrag > 

which implies 

Tdrag ' ^ FBar 

{dE / (it)Bar Rdrag 

Figure 5 shows Tdrag versus a. In computing (Qabs); we have assumed spherical grains, using Mie 
theory and the dielectric functions for silicates from Li & Draine (2001). The specific energy density 
u\ for the ISRF is given in Paper I. From Figure 5 we see that TBar/^drag ‘C 1 for all grain sizes of 
interest when cj ~ wt- Thus, Pdrag • ^ may be neglected in the evolution of q. 


(59) 

(60) 
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3.4. H 2 Formation 


The H 2 formation torque Fhj has a fixed direction in grain body coordinates. When q < I1/I2, 


and 

When q > h/h, 

and 


rH2(^,±) = ± (rH2 • ai) 


h - hq 

h-h 


-,1/2 


TT 


2F{TTl2,k‘^) 


rH2 • ±) = ±|rH2l• 


rH2(Q', ±) = ± (rH2 • %) 


h{q - 1 ) 
h-h 


1/2 


TT 


2F(7r/2,A:-2) 


rH 2 • ±) = ±|rH 2 l• 


(61) 


(62) 

(63) 

(64) 


We estimate the magnitude of the H 2 formation torque as (see Paper II) 

^112 ~ .^lt^H2/'hlrag, gas ; 


(65) 


with 


a;H 2 ~ 5.2 x 10^s“V 


O.l^m Y 
a ) 



En2 
0.2 eV 



nn 


( 66 ) 


/ is the fraction of the H atoms that collide with the grain and depart as H 2 , P is the surface area 
per H 2 formation site on the grain surface, Fhj is the kinetic energy of the departing H 2 molecules, 
and n(H)/nH is the fraction of the gas-phase H that is in atomic form. For a real grain, Fhj varies 
in magnitude and direction on a timescale characteristic of changes in the surface sites. As in Paper 
II, we will ignore this complication and assume that Fhj does not evolve with time. 


Equations (62), (64), and (65) yield 


rH2-^ T-Rar i^H2 

- rsj- 

(ydE / dt)'Qgj^ Tdragigas ^ 


(67) 


Figure 5 shows that rua = ^drag, gasf^/i^T ^ '^Bar when uj ~ cut, so Fhj • ^ may be neglected in the 
evolution of q. 


3.5. Paramagnetic Dissipation 


For steady grain rotation, the torque due to paramagnetic dissipation is given by (Davis & 
Greenstein 1951) 

\"V 


DG = 


UJ 


(ljxB) xB 


( 68 ) 



- 19 - 


Lazarian & Draine (2000) pointed out that is modified by the grain rotation itself. They showed 
that the dissipation rate for a sample rotating in a weak static field B is equal to the dissipation 
rate for a static sample subject to the Barnett equivalent field plus the rotating field B. This 
yields 


X 


// 


_ XO ‘^^2 _ 

1 + 7|B2 sin2 ^TiTs 


~ X0^T2 


(69) 


the second term in the denominator is negligible for the magnetic field strengths encountered in the 
ISM. Since X 0 T 2 for electron paramagnetism substantially exceeds that for nuclear paramagnetism, 
the nuclear contribution may be neglected. Note that equation (69) does not have the same form as 
equation (37) for x” in the case of Barnett dissipation. The rotation of the interstellar magnetic field 
(in the grain’s frame) is in the correct sense so as to be in resonance with the Barnett equivalent 
field lying along ai, so the term involving is zero. The rotating component of the Barnett 
equivalent field, on the other hand, rotates in the wrong direction for resonance. 


The situation is much more complicated in the case of non-steady grain rotation, due to the 
increased complexity of the external field as viewed in the grain’s frame. For example, for an ax- 
isymmetric grain in non-steady rotation, the external field consists of five separate components: a 
constant component along ai, a component that oscillates along ai, and three components that ro¬ 
tate (with different angular speeds) about ai. Even the simplest analysis, in which the components 
of the external field along ai are neglected in comparison with the Barnett equivalent field coiai/xg, 
introduces a complicated dependence on q, since the field strengths and angular frequencies of the 
three rotating components all depend on 


Since we are primarily concerned with radiative torques in this paper, we will ignore the 
complicated g-dependence of the Davis-Greenstein torque, and simply treat it as if the grain were 
rotating uniformly about ai. In this case. 


J 


with 


'Tdg 


2aipa^ 6 

tdg = rr. ^9 ~ 1-5 X 10 yr 


Tdg = + Jsin,^^ 


5X0^252 


aip 


3gcm“3y Y0.I//my \I5Ky \ B 


IjL] 


(70) 


(71) 


Since Tdg • ^ ~ x"^^B‘^-, ^dg ' is suppressed with respect to {dE/dt)Bar by a factor ~ 
(B7g/t(;)2(l -|-(j2TiT 2), rendering it negligible. 


should also be noted that the susceptibilities are not necessarily the same as for the case of a single rotating 
field component. The modified Bloch equations should be solved with all of the rotating field components included. 
It seems unlikely that the susceptibilities would be significantly affected when all of the rotating fields are as weak as 
the interstellar field. However, there could be cases in which the presence of the rotating components of the Barnett 
equivalent field affects the susceptibilities for the components of the interstellar field. 
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3.6. Radiative Torques 


The radiative torque on a stationary grain exposed to a unidirectional radiation field with 
energy density u can be expressed as (Paper I) 




Trad = vra u— (Qr( 0 , ^, /?)) 

ZTT 


where a tilde denotes spectral averaging: 


u = I u\ d\ , 


(72) 


(73) 


and 



u 


A u\ d\ 


Qr(0,d>,/3) = 


\u 


Qr(0,d>,/3,A) A ux d\ 


(74) 

(75) 


We approximate the interstellar radiation field by an isotropic component with energy density 
(1 — 7rad)'Wrad plus a Unidirectional component with energy density yrad^trad, propagating at angle 
■0 with respect to B. We will assume that the starlight anisotropy direction is in the xb — zb plane. 


3.6.1. Unidirectional Component 


We adopt the following transformation between “alignment coordinates” and “scattering co¬ 
ordinates” for the unidirectional component of the radiation field: 


ei = — sin 0 XB + cos 0 zb , 

02 = cos 0 Xb + sin 0 zb , 
ea = YB ; 

see Figure 1. 

The unidirectional component of the radiative torque is then given by 
TracKC, 0, 0, q, ±) = ^AradWadO^A F(^, 0, 0, q, ±)| -b G{C, 0, 0, 9, ±)0 + R(^, 0, 0, Q, ±) J 


where 


= —Qp • ei (sm 0 cos^cos 0 + cos 0 sin,f) 

+Qr ■ (cos 0 cos ^ cos 0 — sin 0 sin .^) + Qp • 63 cos ^ sin 0 , 

G(^, 0,0, g, ±) = Qp • ei sin 0 sin 0 — Qp • 62 cos 0 sin 0 + Qp • 03 cos 0 , 


(76) 

(77) 

(78) 


(79) 

(80) 
(81) 
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and 


,' 0 , 9 , ±) = Qr • ei (cos' 0 cos^ — sin 0 sin,^cos 0 ) 

+Qr • 62 (sin 0 cos ^ + cos 0 sin ^ cos 0) + Qp • 63 sin ^ sin 0 . (82) 


Note that equations (80) through (82) are identical to equations (21) through (23) in Paper II, 
except that here the efficiency factor Qp is averaged over the torque-free motion described in §2.5, 
whereas in Paper II the efficiency factors are averaged over rotation about ai. In order to evaluate 
Qr(^,</>, 0 , 9 , ±), we need to know the values of that correspond to the Eulerian angles 

(a, 0 , 7 ) as the grain executes its torque-free motion. For this, we use equations (7) through (9); 
expressions for the dot products appearing in these equations are given in Appendix B. 

For unidirectional radiation, 


rrad • ^^(0,0,0,9,±) 


r)Trad^rad® A 77(0, 0, 0, 9 , i) 


(83) 


Thus, 

Trad ■ ^ "TSar 

- r\j - 

{dE / (it)Bar '7'rad 

where 


(84) 


Trad = 250 yr 


3 gem 


-3 


1/2 


T \ 


1/2 


100 K 


y/2 / 10-3 X 


O.lfim) V7rad77 


1 I ’'XkJXti- 1 

J V ^*rad J 


(")( 5 ) 


(85) 

We expect 7 rad-f^ ~ 10“^ (Paper II). For the ISRF, A = 1.2/im. Figure 5 shows that TBar/T-ad 1 
when to ~ lot, so Fj-ad • ^ may be neglected in evolving q. 


3.6.2. Isotropic Component 


The torque on a grain exposed to the isotropic component of the radiation field is fixed in 
grain body coordinates and is given by 

nad,0 = ^(l-7rad)Urada2AQr , (86) 

where 

Q|?°.ai = — / dcosQ d0Qp(0,O,/3)-a* . (87) 

dvr Jo 


As with the H 2 formation torque, the components perpendicular to J average to zero during 
the grain rotation. Thus, 


f)l°d(g,±)=±(r)l%.ai) 


h - hq 


-|l /2 


TT 


2F(7r/2,A:2) 


h-h 


( 88 ) 
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when q < I 1 /I 2 and 


rfa°d(^,±) 


±(r 


iso 

rad,0 


•as) 


' h{q - 1) 1 TT j 

/1-/3J 2F(7r/2,A:-2) 


(89) 


when q > Ii/h- 
Also, 

riTd • ±) ~ - 7rad)«rada^A(55?°^ . (90) 

Since (1 — 7rad)Qr° generally less than 7 rad|-f^|, we can neglect • u. 


3. 6 .3. Evaluation of Radiative Torques 

We use the discrete dipole approximation^^ to compute Qr(0,<h,/3) (see Paper I). In this ex¬ 
ploratory paper, we consider only shape 1 from Paper I and a = 0.2/im. We adopt the dielectric 
functions for astronomical silicate from Li & Draine (2001). In order to further reduce the compu¬ 
tational demands, we assume a monochromatic spectrum, with A = 1.2^m (A for the ISRF). In a 
future paper, we will consider a range of grain shapes and sizes and we will fully integrate over the 
ISRF. 

Computations are only performed for $ = 0, since 

Qr(0, 13) = Qr(0, 0 , /3)-eiei-FQr(0, 0 , / 3 )-e 2 (e 2 cos ^+63 sin <I>)-FQr(0, 0 , / 3 )-e 3 (e 3 cos <I>-e 2 sin 

(91) 

We ran the code for 33 values of 0 between 0 and vr (evenly spaced in cos 0) and 32 values of f3 
evenly spaced between 0 and 27r. 

In Figure 6 , we plot F, G, and H versus q for the case that if = 70°, ^ = 30°, and f = 160°. 
For isotropic radiation, Qp° • ai = 1.5 x 10“®, Qp° • a 2 = 1.2 x 10“®, and Qp° • as = —2.2 x 10“'^. 


4. Thermal Averages over q 

Since the Barnett dissipation timescale rear is much shorter than any other relevant timescales 
(except for the timescale characterizing the torque-free motion), we do not follow the evolution of 
q. Instead, we average over a thermal distribution of ( 7 -values when evaluating the torques. We 
consider the following two cases (see §6 for the motivation of these choices). In the first case, 
we assume that the grain is in the positive flip state with respect to ai when q < / 1 // 2 . When 
q > / 1 // 2 , we average over both flip states (with respect to a 2 when q = I 1 /I 2 and with respect 


use the code DDSCAT.SalO, which is available at http://astro.Princeton.EDU/~draine/. 
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to as when q > Ii/ 12 ). A quantity A averaged in this way is denoted by (A)_|_. (We assume that 
A has already been averaged over torque-free motion and do not include a bar to indicate this 
averaging.) The second case is identical to the first, except that the grain is in the negative flip 
state with respect to ai when q < h/h] such averaging is denoted by (A)_. 


4.1. Density of States 


In order to take thermal averages over q, we need to evaluate the density of states in q (for a 
grain with constant J). To this end, we examine the grain’s trajectories in phase space. In general, 
the phase space for a rotating grain has six dimensions: one for each Euler angle plus one for each 
conjugate momentum. However, when the angular momentum is specified this reduces to just one 
angle plus its momentum. 

The Lagrangian L is just the kinetic energy: 




t9 9 t9 ■ 9 ■ 9 t9 9 • 9 

2 J COS 7 j^sin asm 7 J^cos asm 7 


2/1 ^ 2/2 ^ 2/3 


(see eqs. 29 through 31). The iVi are related to the Eulerian angles as follows: 

uoi = ( cos 7 -|- d , 

<^2 = C sin a sin 7 -|- 7 cos a , 

and 

UJ 3 = C, cos a sin 7 — 7 sin a . 

Using these relations, we rewrite the Langrangian as 

r [t , u T \ 2 1 _ hhh^^ _ 

“ 2 hh ^ 2 [{h - h)h + h{h - h) cos 2 a] • 


(92) 

(93) 

(94) 

(95) 

(96) 


The conjugate momentum is given by 


dL — I 1 I 2 I 3 A 

da [{h - l 2 )h + h{h - h)cos‘^ a] 


(97) 


Expressing the second term in equation (96) for the Lagrangian in terms of pa and setting L equal 
to the energy E, we find 


Pa = J 


h {h - l2q)+ h {h - h ) cos^ 


a 


1/2 


h{h-h) + h{h-h)cos^a 


( 98 ) 


Sample phase space trajectories are plotted in Eigure 7 for the case that Ii ■. I 2 : = 2 > : 2 : 1 . 
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The number of states with energy between Ei and E 2 is proportional to the area in phase 
space bounded by the trajectories for Ei and £^2- Thus, the number of states with g-value between 
1 and q is proportional to the following quantity: 



- hq) + h{h - h) cos^ 


a 


h{h-h)+h{h-h)cos^oc 


1/2 


where 

, q < h/h 
, q > h/h 

Note that s ranges from 0 to 1 as ranges from 1 to h/h- 


ai — 
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1 —1 

\h(hq-hA 

cos 

h{h—h) 
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Since the density of states in the variable s is constant, we will evaluate thermal averages by 
integrating over s. Figure 8 shows g as a function of s for a grain with Ii : I2 ^ -^3 = 3 ; 2 ; 1 and 
for shape 1 from Paper I. 


4.2. Thermal Averages for the Torques 


The average of a function A{q, ±) over a thermal distribution of q at constant J is given by 

fo" dsA[q{s),±] e^p[-q{s)h/2IikTd] + ds0.5 {A[g(s),+] + A[q{s), -]} e^p[-q{s)h/2hkTd] 


^ J ds exp[—q{s)h/2IikTd] 

where the value of s when q = h/h is given by 


Sr = 1-sin 


-1 


TT 


/l(/2-/3)l'/" 


h{h-h) 
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If A{q) does not depend on the flip state, then we will simply denote the thermal average over q 
by (^). 


The Barnett, drag, and Davis-Greenstein torques are all invariant under a change in flip 
state. With our simple approximations, the Davis-Greenstein torque does not even depend on q, 
so (Tdg) is given by equation (70). We also have approximated the Barnett and drag torques 
as simply proportional to q. Thus, (Tb) is given by equation (50) with q replaced by {q) and 
(rdrag) = -(9)J/'rdrag- In Figure 9 we plot {q) versus h/{2IikTd). As h/{2IikTd) 00, {q) 1 

and as J‘^/{2IikTd) —> 0, {q) fg dsq(s). 


The H2 formation torque does depend on flip state and its thermal average is given by 


where Ci = {C) and 


(rH2)± = ± 


h^n2 

"^dragjgas 



ai) CiJ 
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h-hq 
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TT 
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, q < h/h 
, q > h/h 


(103) 


(104) 
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Note that there is no contribution to (rH2)± from q > I 1 /I 2 , since we average over flip states 
for these q, with perfect cancellation. In Figure 9, we plot Ci versus /{2IikTd); Ci ^ 1 as 
J^I{2hkTd) ^ 00. 

The thermal average of the torque due to the isotropic component of the radiation field is 
similar to that for the H2 formation torque: 

ra)± = ±^(l-7radKada"AQr-aiCiJ . (105) 

The thermal average of the torque due to the unidirectional component of the radiation field 
requires the thermal averages {F)±, {G)±, and {H)±. It is clear from Figure 6 that these do not go 
to zero when averaging over flip states. (It is not always the case, however, that H{q) for the two 
flip states are nearly identical.) In Figure 10 we plot {F)±, {G)±, and {H)± versus /{2IikTd) for 
the case that V’ = 70°, ^ = 30°, and (p = 160°. 


5. Equations of Motion 


Suppose that during some time interval At [during which the grain has angular momentum J 
defined by (J, p)], if the grain is in a flip state with respect to ai, then the probability that it is 

in the positive flip state is /+, and the probability that it is in the negative flip state is /_ = 1 — /+. 
In this case, the thermal average of a quantity A is given by 


{A) = {A)^U + {A)_f_ . 

Prescriptions for /+ and /_ are developed in §6. 

The following equations of motion result from equation (11): 

# / , 7rad^^rada^A(G)(^,()),'0, J) 

— = (g)IlB +- r,T ■ c - 

at 2Jsm4 

^ ^ 7 rad^rada^A(F)(^, (/>, j;, J) _ siu ^ COS g 

dt 


and 


2J 


I, r)7rad^radQ' A(77)(^, (/>,'0, T) 

dt 2 


'TDG 


J 


"^drag 


+ 




1 


+ ^(1 - 7rad)^irada^AQr™ 
L^^dragjgas ^ 


aiCi(/+-/_)- 


J sin^ 
tdg 


(106) 


(107) 

(108) 


(109) 


These equations are similar to equations (24) through (26) in Paper II (except note the sign error 
involving Hb in eqs. 10 and 24 in Paper II). The dependence of (F), (G), and {H) on J results 
from the presence of J in equation (101). 
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In Figure 11 we plot the dynamical times associated with the terms involving (-F)+, (G)+, and 
{tf, tg, and th, defined in the figure caption), as well as the magnetic precession timescale 
the Barnett dissipation timescale rear) and the thermal flipping timescale rtf (see §6), 
versus for the case that = 15K, ip = 70°, ^ = 30°, and (p = 160°. As mentioned in 

§4, TBar is much shorter than the other timescales (although note that Tp, tq-, and th are shorter 
when the entire ISRF is adopted). 

Since ((g)llB)~^ is much shorter than Tp-, tq, and th for all J‘^/{2IikTd) in Figure 11 (and 
the radiative torques are generally expected to be more important than the Davis-Greenstein and 
H 2 formation torques; see Paper II), we expect that we will generally be able to average the 
equations of motion over grain precession. In this case, only equations (108) and (109) remain, and 
(F)(^, (p, Ip, J) and (p, ip, J) are averaged over the motion in (p. Thus, we define 

{F)i{^,iP,J) = - r dcp ^ {F)±{i,cP,iP,J) (110) 

T/) Jo dt 

and ^ 

r dcp ^ {H)±i^,cP,iP,J) , (111) 

'^(j) Jo 

where 

d(p\d(p / dt\~^ . (112) 

When Tg dcp/dt is nearly independent of (p and (T)± and (i7)± have the form adopted 

in Paper II (these were denoted F and H in Paper II). If \tg\ < (((?)Hb)~^ and tg{(P) undergoes 
sign changes, then the grain might not even undergo complete cycles of 0 to 27r in (p. However, we 
expect that this will only occur for ^ ^ 0 or tt. 

When averaging over the motion in <p, we replace {F){^, (p, ip, J) and (p, ip, J) in equations 

(108) and (109) with {F)i{^,iP, J)U + {F)t{^,iP, J)f. and {H)^4^,iP, J)U + {H)t{^,iP, J)f_, 
where the evaluation over the thermal averaging cases is done after the averaging over the motion 
in (p. When rtf < ((g)!!^)”^, these evaluations should actually be done in the reverse order, but in 
that case it would not be possible to construct an interpolation table for (T)± and (77)^ prior to 
running the evolution code (see §6). Since it is almost always the case that tg {{q)^B)~^, this 
reversal of order will not introduce significant error. 

The behavior of {F)'^(^,ip, J) in the limit ^ ^ 0 or tt is somewhat subtle; see Appendix C for 
a detailed discussion. 


6. Algorithm for Evolving Dynamics With Thermal Flips as the Only Stochastic 

Element 

As mentioned in §1, we expect stochasticity due to random gas atom impacts and H 2 for¬ 
mation to be less important than random thermal flipping due to Barnett relaxation. Thus, as 
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a first approximation, we will evolve the grain dynamics without taking the former stochasticity 
into account. In a future paper we will describe a more detailed analysis that accounts for this 
stochasticity in an approximate way. 


In order to ensure a consistent treatment of the stochastic flipping throughout the entire 
dynamical evolution, we integrate the equations of motion using a constant time step size At. 
We do not expect that there would be any significant biases if a scheme with a variable step size 
were employed. However, when we generalize further (in a future paper) to include stochasticity 
associated with gas atom impacts and H 2 formation, we will adopt Langevin equations, for which 
constant step sizes are preferred. Even with a constant step size, the computer time spent running 
the evolution code is insignificant compared with the time devoted to the scattering calculations 
and construction of interpolation tables for (E)^ and (77)^. 


For sufficiently large J, the grain will rarely have q > Ii/l 2 - In this limit, we can consider the 
grain to always be in one of the flip states with respect to ai, with some probability per unit time 
that a flip will occur. We adopt the simple approximation suggested by Lazarian & Draine 
(1999a); 
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During time interval At, the grain spends some fraction /o of the time in its original flip state, 
and a fraction (1 — /o) in the opposite flip state. Ideally, we would construct the distribution 
function for /q and randomly select an /q from this distribution for each time step. However, it 
is not clear how to construct this distribution function. The expectation value /same(At) for the 
fraction of the time step that the grain will be in its original flip state is easily obtained; 


/same (At) 


T r dt = l\i + «>^P(-A«/y)sinh(A(/n,) 

At Jo ^ 2 [ At/rtf 


(115) 


where 

Psame(At) = exp(-At/rtf) cosh(At/rtf) (116) 

is the probability that the grain is in the same flip state after a time At (see Appendix D for the 
derivation of this result). As a simple approximation, we could set /o = /same for each time step. 
This would represent “typical” behavior”, but would deny the grain the chance of ever spending 
a significant period of time in one particular flip state. Such a “run” in one flip state could occur 
in reality, and could have the important consequence of allowing a grain that would otherwise be 
stuck in the thermal rotation regime to escape to suprathermal rotation. 
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Thus, we adopt the following slightly more complicated prescription. The probability that the 
grain undergoes no flips during the time step is 


Po = exp(-At/rtf) 


(117) 


We take the grain to remain in its original flip state throughout the time step with probability Pq. 
If the grain does not remain in its original flip state throughout the time step, then we assume 
that the grain spends a fraction fg of the step in its original flip state, where fg is determined from 

-Pq + (1 “ Po)fs = /same- 


fs = 


/si 


-Pn 


(118) 


1-Po 

Thus, if the grain is in the positive flip state at the start of the time step, then /_|_ = 1 and /_ = 0 
with probability Pq; otherwise, /+ = fg and f- = 1 — fg. 


The flip state at the start of each time step is a stochastic variable. As mentioned above, we 
take the grain to remain in its original flip state throughout the entire step with probability Pq. In 
this case, the flip state at the start of the next step is taken to be the same as it is in the current 
step. Otherwise, the flip state at the start of the next step is taken to be the opposite of the flip 
state at the start of the current step with probability Pf, determined from 1 — Psame = (1 — Po)Pf: 

Pf = . (119) 

1 - Po 

Table 1 summarizes the 3 branches which can be taken on each time step. 

In the limit of low J, the grain does spend a significant fraction of the time with q > / 1 // 2 . 
Furthermore, the grain undergoes rapid flipping in this limit, so that the grain spends half of its 
time with q > I 1 /I 2 in the positive flip state with respect to as and and half of the time in the 
negative flip state. This motivates the definition of our two thermal averaging cases, given at the 
beginning of §4. 

If At > r 0 , then we average over motion in (/>; otherwise we do not. If (/>-averaging mode 
applied in the previous time step but does not apply in the current step, then a value of 4> is selected 
randomly at the start of the current time step. Since we expect that i^-averaging mode will always 
apply, we tabulate {F)±[cos J‘^/{2IikTd)], {H)±[cos^, /{2IikTd)\i and r 0 [cos J^/(2/iA:rrf)] for 

a given tp before running the dynamical evolution code and interpolate as needed. If non-(/>-averaged 
quantities are needed at any point, then they are calculated within the evolution code. 


7. Results 

As an example of our new method for evolving the grain dynamics including the effects of 
thermal fluctuations and thermal flipping, we consider shape 1 , a = 0 . 2 ;um, a monochromatic 
radiation field with A = Aisrf = 1 . 2 /im, Urad = ^^rad.iSRF) 7rad = 0.1, T = 100 K, nn = 30cm“^, 
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= 15 K, i? = hfiG, p = 3g cm“^, and ip = 70°. Because we adopt a monochromatic radiation field 
with A = 1.2/im (as opposed to using a realistic spectrum for interstellar starlight), this example will 
be illustrative of a case where the radiative torques—while still dynamically very important—are 
relatively weak (corresponding to a ~ O.l/um grains with a realistic starlight spectrum). 

For simplicity, we ignore the torques due to infrared emission (which for this case is negligible 
compared to gas drag) and H 2 formation. In future work we will ascertain the importance of the 
H 2 formation torque. 

In Figure 12, we plot the trajectory map calculated using the method of Paper II, in which 
it is assumed that J always lies along ai (parallel or anti-parallel). The dotted (solid) trajectories 
originate at J/Iilvt = 30 (0.5) and the upper (lower) half-plane is for the positive (negative) flip 
state. 

If we had included this map in Paper II, then it would have looked quite different. In Paper 
II, ^ was defined as the polar angle of ai in alignment coordinates, whereas here we have defined ^ 
as the polar angle of J in alignment coordinates. Henceforward, we will write for the angle as 
defined in Paper II. When the grain is in the positive flip state with respect to ai, but in 

the negative flip state, Cii = In Paper II, negative flip states were indicated by assigning a 

negative value to J; here J = |J|, which is always positive. Since trajectory maps in Paper II were 
drawn in terms of cosc^ip the lower half of the map in this paper (i.e., below the line J = 0) must 
be reflected in the line cos.^ = 0 in order to look like the maps in Paper II (and vice versa). 

Note the near symmetry of the map under reflection in the line J = 0. The maps in Paper 
II also appear somewhat symmetric, though to a lesser degree than here. (In the maps in Paper 
II, the symmetry is with respect to a 180° rotation.) Evidently the radiative torque changes less 
under a grain flip for a monochromatic radiation spectrum with A = 1.2|Um than for the ISRF (at 
least for shape I). 

From Figure 12, we see that there are two stable stationary points (“attractors”), with cos^ = 1 
and J/Ii^T = 17.7 (14.8) for the positive (negative) flip state. There are also two unstable 
stationary points (“repellors”), with (J/Iiwt, cos ^) = (8.3,0.76) for the positive flip state and 
(J/ZitJT, cos.^) = (4.2,0.74) for the negative flip state. 

Crossover points, where the trajectories cross J = 0, are located at cos .^11 = —1, -0.75, 0.77, 
and 1. (Those at cos.^n = and 1 are “crossover repellors”. Only one trajectory passes through 
each of these crossovers; thus, they are not physically relevant.) Suppose J really were constrained 
to lie along ai, for all J. In this case, a grain could enter the crossover at cos^n = —0.75 from 
above, i.e., in the positive flip state. It would then emerge in the negative flip state, still with 
cosc^n = —0.75, in which case cos^ = 0.75. The direction of ai remains unchanged during the grain 
flip, while J flips. In other words, the grain’s orientation does not change; it stops spinning, and 
then begins spinning in the opposite sense. Some of the trajectories emerging from this crossover 
proceed to the crossover at cos .^11 = 0.77 (cos^ = —0.77) and others lead to the attractor at 
cosc^n = — 1 (cos^ = 1). Trajectories emerge from the crossover at cos.^n = 0.77 with the grain 
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in the positive flip state, with some heading towards the crossover at cos^ii = —0.75 and others 
heading towards the stationary point at cos^ii = 1. Thus, in the terminology of Paper II, this is a 
“semi-cyclic” map. 

The dynamical picture changes in two distinct ways when the formalism developed in this 
paper is adopted. First, since we no longer enforce ai || J, the structure of the torques as functions 
of J, and flip state are different. Thus, even in the absence of thermal flips, we can expect 
qualitative changes in the structure of the trajectory map at low J. Second, it is no longer possible 
to construct strictly deterministic trajectory maps, because of the random thermal flipping. In 
order to gauge the importance of these two effects, we will progress in stages, first suppressing the 
thermal flips. 

In Figure 13, we display (F)!^(^) and (i7)+(^) for J ^ oo, for J ^ 0, and for an intermediate 
T; {F)t_{0 and {H)t_{0 are nearly identical to (F)^(^) and (i7)+(C) at low J. The low-J versions 
of these functions are quite different from the high-J versions. Since the locations of the stationary 
points and crossovers depend on the structure of these functions, this implies that the locations of 
crossovers and low-J stationary points could differ from the locations inferred using the analysis of 
Paper IF Indeed, there could be crossovers or stationary points not even indicated by the analysis 
of Paper II. Alternatively, such points found using the Paper II formalism might turn out not to 
exist after all. In Paper II, we showed that crossovers occur when (F’)'^(^) = 0. Since {F)'^ only 
has zeros at cos^ = ±1 as J —> 0, we may expect to no longer find crossovers at cos^ii = —0.75 
and 0.77. 

Figure 14 is the trajectory map (with J/Iiwt = 30 initially) constructed using the method 
developed in this paper, except that thermal flipping is prohibited: a grain which begins the 
trajectory in the positive flip state has /+ = 1 and /_ = 0 throughout (and vice versa for the 
negative flip state). Since thermal flipping is not included, these are “deterministic” trajectories: 
they are fully determined once the initial conditions (including flip state) are specified. Note the 
dramatic change in the map: two new attractors have appeared, at ( J/Ii^t, cos ^) = (1.40, —0.57) 
(positive flip state) and (0.71,-0.49) (negative flip state), and the crossovers have disappeared. 
The map is now non-cyclic. 

The low-J attractors arise because (JI)^ are smaller in the low-J limit than in the high-J limit 
due to orientational averaging—at small J, thermal fluctuations allow the grain to spend a large 
fraction of the time with large angles 7 between J and ai (see Fig. 13; curves labelled “b” are for 
the value of J at which the new attractor appears for the positive flip state). We have found that 
the low-J attractors remain even if the radiation field intensity is increased by a factor 100. Since 
the low-J attractors appear even in the absence of thermal flipping, they clearly are not a thermal 
trapping effect (Lazarian & Draine 1999a; see §1). The radiative torque simply drives the grain to 
a state with low J. 

When a grain undergoes a thermal flip, J remains fixed in space (i.e., J and ^ are unchanged) 
and ai flips. In Figure 12 or 14, the point representing the grain’s state is reflected in the line J = 0 
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when the grain thermally flips. Thus, thermal flips do not carry grains across radiative torque- 
induced crossovers. Because of the near symmetry of the trajectory map, if a grain undergoes a 
thermal flip while approaching one of the low-J attractors, it will in most cases land on a trajectory 
approaching the other low-J attractor. Similarly, a grain that thermally flips while approaching one 
of the high-J attractors will usually land on a trajectory approaching the other high-J attractor. 
Due to the presence of the repellors at cos^ ~ 0.75, there are some cases where a grain that is 
approaching one of the low-J attractors can be redirected towards one of the high-J attractors (or 
vice versa), but cos^ must start off quite close to 0.75 for this to happen. 

Thus, the description of the grain dynamics using deterministic trajectories is an excellent 
approximation in this case. This is not necessarily true generically, but relies on the high degree of 
symmetry in the map. Thermal flipping does, however, change the appearance of the trajectories. 
For sufficiently small J, the grain undergoes rapid flipping, so the trajectory map must be symmetric 
in reflection about the line J = 0. In Figure 15, we display trajectories constructed using our full, 
stochastic evolution algorithm (for J/Ji^t = 30 initially). For trajectories in the upper (lower) 
half-plane, the grain is initially in the positive (negative) flip state. The expectation values /+ 
and /_ are not explicitly indicated at any point along the trjactories, but it is clear that /+ ~ 1 
when J/Iicut > 5 (for trajectories that start in the positive flip state) and f+ ^ f- ^ 0.5 when 
J/IicuT < 5. The transition from virtually no flipping to rapid flipping occurs when J/Iicut ~ 5, 
as expected from Figure 11. We have displayed trajectories resulting from one iteration of our 
stochastic evolution code; other realizations look nearly identical. 

The low-J attractors fall in the rapid-flipping regime; the rapid flipping shifts their position to 
{J/Iiujt,cos^) = (1.14,-0.50). Although Figure 15 gives the impression that there are two low-J 
attractors, these are really the same, because of the rapid flipping. 

In Figure 16, we plot the map for trajectories starting with J/Iicut = 1. The grain is always 
in the rapid-flipping regime for the low values of J in Figure 16. All of these trajectories land on 
the new, low-J attractor (except for physically irrelevant trajectories with cos^ = 1 initially). We 
have found that when J/Iicut < 2 initially, all trajectories land on the low-J crossover. Thus, if 
a population of real interstellar grains were characterized by a trajectory map similar to this one, 
then the high-J attractors might not be relevant, since the grains are unlikely to be in an initial 
state of suprathermal rotation (although only a very mild initial suprathermality is needed for some 
grains to end up at the suprathermal attractors). 

We have also analyzed the grain dynamics when ^ = 0 and 30°. The method of Paper II yields 
a semi-cyclic map with two suprathermal attractors when ■0 = 0 and a cyclic map when 0 = 30°. 
In both of these cases, the method developed in this paper yields a non-cyclic map with a low-J 
attractor. 
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8. Summary 

In Paper II, we derived equations of motion for a grain exposed to radiative torques, under 
the simplifying assumption that the angular momentum J is always parallel to the principal axis of 
greatest moment of inertia ai. We found that radiative torques can spin grains up to suprathermal 
rotational speeds and (in combination with the Barnett torque) align the grains with the interstellar 
magnetic field. In addition, radiative torques can drive grains into periods of thermal rotation, when 
they might undergo a type of crossover in which the direction of J reverses. 

Because of internal fluctuations, ai and J are not necessarily parallel when J is small. Thus, 
we were unable in Paper II to follow the grain dynamics during periods of thermal rotation. Con¬ 
sequently, we were unable to ascertain the role of crossovers in the dynamics. 

Here, we have generalized the treatment of Paper II, relaxing the constraint that ai || J and 
including the effects of thermal fluctuations. Our principal results are: 

1. We provide a detailed description of torque-free rotation (§2.5), considering the most general 
grain shape, with no degeneracy in the eigenvalues of the inertia tensor. The rotation is completely 
characterized by the angular momentum J, the rotational energy E, and the “flip state” (§2.5.2). 
For convenience, we define a dimensionless measure of the energy at constant J, q = 2IiE/j‘^. 

2. We derive approximate expressions for the external torques, averaged over torque-free rotation 
(which occurs on a much shorter timescale than the dynamical timescales associated with the 
torques), as functions of J, q, and flip state. We consider the torques due to the Barnett magnetic 
moment, gas and IR emission drag, H 2 formation, paramagnetic dissipation, and radiation. 

3. We make an order-of-magnitude estimate of the Barnett dissipation rate (§2.5.4) and find that 
internal dissipation dominates the external torques in determining the evolution of q. 

4. We derive the density of states in the parameter q and find expressions for the thermal averages 
of the torques for various cases of interest (§4). 

5. We develop an algorithm for evolving the grain dynamics (§6). The algorithm allows for random 
thermal flips. 

6 . As an illustration of the method, we calculate a trajectory map for one particular case (§7). 
Using the method of Paper II, we find two attractors (characterized by suprathermal rotation) and 
two crossover points. The more complete treatment developed here yields no crossover points and 
a third attractor, characterized by thermal rotation and rapid flipping. Thus, there are qualitative 
differences introduced by the proper inclusion of thermal fluctuations and flipping: whereas the 
Paper II analysis produces a semi-cyclic map, the new analysis produces a non-cyclic map. 

In future work, we will investigate more cases, in an effort to determine whether or not non- 
cyclic maps with low-J attractors are a generic feature when the analysis developed here is applied. 
For this, it will be crucial to employ the full ISRF, rather than a monochromatic spectrum with 
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A = 1.2^m. In addition, we will explore the effect of stochastic gas atom impacts and H 2 formation 
events on the dynamics. Ultimately, we will conduct a systematic study of grain alignment by 
radiative torques, considering a range of shapes and sizes. 


This research was supported by NSERC of Canada, NSF grant AST 9988126, and an NSF 
International Research Fellowship to J. C. W. We are grateful to A. Lazarian for many stimulating 
discussions. 


A. Angles in Scattering Coordinates 

When ai • 63 = 0, $ = 0 if ai • 62 > 0 and <I> = tt if ai • 62 < 0. 
When sin 0 = 0, 


(3 = 2 tan 


/I - a 2 • 62 cos 0 


(Al) 


a2 • 03 

If a2 • 63 = 0 , then /3 = tt — 0 if a2 • 62 = —1 and /3 = 0 if a2 • 62 = 1. 

When sin 0/0 and a2 • 63 cos <I> = a2 • 62 sin < 1 >, /3 = 0 if a2 • ei < 0 and /3 = tt if a2 • ei > 0 . 


B. Dot Prodncts Needed in §3.6.1 


with 


a* • ei = cos / Zb • a* — sin / xb • aj , 
a* • 62 = cos / XB • a* + sin / zb • aj , 

yB ' 1 

xb • ai = cos ^ cos (j) sin C sin 7 + sin (j) cos C, sin 7 + sin ^ cos cf) cos 7 , 

ZB • ai = cos ^ cos 7 — sin ^ sin C sin 7 , 

Yb • ai = cos ^ sin cf) sin C, sin 7 — cos 4> cos C sin 7 + sin ^ sin (j) cos 7 , 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B 6 ) 


XB-a 2 = cos^ cos 4 * (cos a cos C — sin a sin (4 cos 7 )—sin(^ (cos a sin + sin a cos C cos 7 )+sin,f cos (/isina sin 7 

(B7) 

Zb • a 2 = — sin ^ (cos a cos — sin a sin (4 cos 7 ) + cos ^ sin a sin 7 , (B 8 ) 

and 

Yb •^2 = cos ^ sin (j) (cos a cos C — sin a sin Q cos 7 )+cos (p (cos a sin + sin a cos C cos 7 )+sin sin p sin a sin 7 

(B9) 
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C. Calculation of {F)‘^{^) when sin^ 1 


When sin C « 1, (Qr)±(?)</>) ~ (Qr)±(^ = 0 or tt) is nearly constant as the grain orbits in 4> 
and, from the definitions of F and G (eqs. 80 and 81), we can write 


and 


where 


Fl^ 


{F)±{i, 4>) = F± cos^cos((/> - (t)o) - sin^ 

(G')±(^,(/>) =-F 2 sin((/)-(/)o) , 

((Qr) 


± • es 


((Qr)± •eisinV’-(Qr)± •e 2 CosV') + 

F[ = (Qr)± •eicosV' +(Qr)± •e 2 sinV' , 

and the angle (po is determined by 

sin 00 = (Qr)±-e3/Fi 

cos00 = (-(Qr)± • ei sm0 + (Qr)± • ^2 cos0)/f£ 
The equation of motion (107) for 0 becomes 


1/2 




1 - 


g 


sin^ 


sin(0 - 0o) 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 

(C7) 


where 


»-^)W^FBx. 0 -(^)(^)(^)(^)’'^(^)-‘ ; (C 8 , 

in equation (C8), we have taken {q) ^ 1, yj-ad = 0.1, Tgas = 100 K, xo = 3.3 x 10“^, and B = 5 fiG. 

A transition from periodic motion in 0 to no motion in 0 occurs when sin^ = g. When 
sin.^ > g, the motion in 0 is periodic with period 


r</ = 


27r 


{q)^b 


“ 

/ \ 2 l 

1 - 

{ ■\) 


\sm^J 


(C9) 


and cos(0—0o) averages to zero over the motion in 0. When sin.^ < g, the dynamics is characterized 
by a stable stationary point 0<j given by 


sm(0^-0o) = sin^/g 

cos{(ps - 0o) = [l - (sin^/c/)^] 


Thus, 


{F)U0 = 


—Fj_ sin^ 


—Fj. sin^ + cos.^ 1 — {sin^/g^ ,s\\i^<g, sin,^ 4/1 


, 5 < sin^ <C 1 


(CIO) 

(Cll) 


(C12) 
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As discussed in § 6 , we construct interpolation tables in cos.^ and /{2IikTd) for and 

(F)^, to be used in evolving the grain dynamics. However, we use equation (C12) to directly 
evaluate (F)^ when ^ lies within one interpolation zone of 0 or tt. Otherwise, it would be impossible 
to reproduce the small-sin,^ structure with a reasonable number of interpolation zones. Thus, we 
also construct a table of /{2IikTd)\. To ensure continuity at = the first interpolation point 
away from ^ = 0 or vr, we estimate Fj_ ~ — (F)^ (?i)/sin6- 


D. Derivation of Fgame 

The probability that exactly n flips occur during At is 

Pn = , (Dl) 

n! 

where k = At/rtf. The probability of zero flips or an even number of flips is 



(D2) 
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E. Glossary of Notation 


symbol 

a 

ai, 3.2, a3 

h,h,h 

J 

q 

Td 

TBar 

A 

{A)±iJ,Td,^A) 

{A)i{J,Td,0 

Po{At) 

Psame(At) 

/same (At) 

/s(At) 

Qr(0,‘b,/3,A) 


significance 

radius of sphere of equal volume, a = ( 3 F/ 47 r)^A 
grain principal axes (eigenvectors of moment of inertia tensor) 
eigenvalues of moment of inertia tensor (/i > I 2 > I 3 ) 
angular momentum of grain 

rotation rate for a sphere in thermal equilibrium with the gas (eq. 10) 

E/(JV2/i) (eq. 13) 

vibrational temperature of grain 

at fixed J, timescale for Barnett damping 

of rotational kinetic energy in excess of J^/2Ii (eq. 43) 

time-average of dynamical variable A over torque-free motion 

with E/(J^/2/i) = q and in ± flip state (eq. 32) 

average over radiation field spectrum (eqs. 74, 75) 

thermal average of dynamical variable A 

for grain with given J and Td, in ± flip state with respect to ai (eq. 101) 

thermal average of variable A for grain with given J and Td, in ± flip 

state with respect to ai, averaged over motion in 4 > (eqs. 110, 111) 

probability per unit time that grain will change its flip state (eq. 113) 

probability that grain undergoes no flips during interval At 

probability that flip state at end of interval 

will be same as flip state at beginning of interval (eq. 116) 

expectation value for fraction of time that grain 

will spend in initial flip state during interval At (eq. 115) 

expectation value for fraction of time spent in original flip state 

for grain which flips at least once during interval At (eq. 118) 

radiative torque efficiency vector for grain orientation (0, <h, (5) 

relative to direction of incident radiation with wavelength A (eq. 72) 
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Table 1. Branching Probabilities 


probability 

fo 

final flip state 

Po 

1 

original (no flips) 

(l-Po)(l-Pf) 

fs 

original (even ^ of flips) 

(1 - Po)Pi 

fs 

opposite (odd ^ of flips) 
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-> 


B 



ang. momentum frame Xj,yj,Zj 
grain axes apiig.ag 
Euler angles a, y 





align, frame Xg,yg, 
scatt. frame e^Cg, 



grain in scattering frame 


Fig. 1.— Orientation of the “angular momentum” (xj,yj,zj) and “scattering” (ej) coordinate 
frames in the “alignment” (xb , ye, zb ) coordinate frame and orientation of the grain axes a* in the 
angular momentum frame and scattering frame. 


<tsi <cu 
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t/Pt 


Fig. 3.— Components of the angular velocity along grain principal axes, for Ii : I 2 ■ I 3 = 3 : 2 : 1 
and q slightly less than h/h] + (-) indicates the positive (negative) flip state with respect to ai. 
(jj 2 does not depend on the flip state. 








0.001 


0.01 


a (/xm) 


0.1 


1 


Fig. 5.— Timescales for Barnett relaxation (rear)) the drag torque (rdrag), the H 2 formation torque 
(thj), and radiative torques (Trad) for silicate grains, assuming uj = wt; -^1 ~ 2 / 3 , = 15K, 

p = 3gcm“^, T = 100K, u-h = 30cm“^, ttrad = "Wisrp, / = 1, ^ = lOA, = 0.2eV, n(H) = u-h, 

and 7 rad-f^ = 10“^. Both thj and Trad oc and, coincidentally, are nearly identical for the above 
canonical conditions. 



Fig. 6.— F, G, and H 
and positive (solid) an( 
Aisrf = 1.2nm was adc 
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Fig. 8.— q versus s for a grain with Ji : /2 • .^3 = 3 ; 2 : 1 and for shape 1 from Paper I. 











-48- 



10-2 0.1 1 10 102 102 104 

JV(2I,kT,) 


Fig. 10.— {F)±, {G)±, and {H)±, normalized to their values when q = 1, for shape 1, one particular 
angular momentum direction (^, (f>) and starlight anisotropy direction ip. 





JV(2I.kTj 


Fig. 11.— The Barnett dissipation timescale TBar, the thermal flipping timescale rtf, and timescales 
associated with various terms in the equations of motion: tf = 2J/(7i.ad'i*rada^^(-^)+)) tg = 
2J sin,f/(7rad'Wrada^A(G)+), and th = 2J/(7i.aditrada^A(F)+), for shape 1, a = 0.2;um, = 15K, 

f/, = 70°, ^ = 30°, and (j) = 160°. 





Fig. 12.— Trajectories computed using the method of Paper II (i.e., constraining J to be parallel 
or anti-parallel to ai) for if) = 70° (see text for specification of other parameter values). Dotted 
(solid) curves originate at J/Iicjt = 30 (0.5). Attractors (repellors) are indicated by open circles 
(filled triangles). 
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-1 - 0.5 0 0.5 1 


cos ^ 

Fig. 13.— {F)^ and (H)'^ versus cos^ for shape 1 and V’ = 70°. a) J ^ oo; b) J/huJT = 1-7 (for 
Tgas = 100K and = 15K); c) J ^ 0. 



Fig. 14.— Trajectory map (with same parameters as for Figure 12) computed using the method 
developed in this paper (i.e., relaxing the constraint that J must be parallel or anti-parallel to ai), 
except that thermal flipping is prohibited (see text for details). 














